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I \ ABSTRACT 

. This study deals with vibroacoustics under heavy fluid loading conditions. Con- 

siderable attention has been and remains focused on this subject not only because 
. industry is very concerned but also because of mathematical difhculties that make 

^ ' the numerical resolution of the problem very difficult. 

1— H , ft was recently observed [8| in a numerical study on a high order perturbation 

^ ' method under heavy fluid loading that a loaded vibrating plate results in a fre- 

^ , quency shift of the in vacuo single resonance (in both the real part because of 

• i-H ' the fluid added mass and the imaginary part because of energy lost by radiation 

into the fluid) as well as increase in the number of the resonance frequencies : as 
' a result of the loading, each single in vacuo resonance frequency of the structure 

CIh. is transformed into a multiple resonance frequency. Here we show that this phe- 

nomenon is said to be an extension to the heavy loading condition of the Sanchez's 
, classical result that have established that in the case of a light loading conditions 

^ ' "the scattering frequencies of a fluid loaded elastic structure (ie the resonance fre- 

, quencies) are nearly the real eigenfrequencies of the elastic body alone and the 

complex scattering frequencies of the fluid with a rigid solid" . 

Using classical results in the framework of the theory of entire functions, it is 

, established that a single resonance of a simply supported fluid loaded rectangular 

, plate is transformed into an infinite number of resonances. 

oo . 
o ■ 

> : 1 INTRODUCTION 

It is well known that the need to accurately describe fluid structure interactions in- 
creases the computational cost of the various numerical methods used. Most numerical 
methods can be used to deal with the resolution of linear systems of simultaneous equa- 
tions. When the loading conditions are taken into account, the size of the matrices 
involved increases considerably, and they become full and frequency dependent. Any 
method of reducing these drawbacks is therefore most welcome; and one of the best 
methods available is that based on the perturbation approach. Asymptotic analysis can 
be conduced when the loading is light, as in the case of a metal structure surrounded by 
air. This approach consists in introducing a small parameter e, which is the ratio between 
the surface mass density and the fluid density. Using perturbation expansions [9], one 
can then construct an approximate solution, based on the in vacuo eigenmodes, which 
is not only easier to calculate than the solution of the exact problem but also leads to a 
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better understanding of the various phenomena involved p]. 

The hght fluid loading perturbation method (involving high density structures in 
contact with a small density fluid such as aluminum plate in contact with air) was recently 
extended [8] to cases where the perturbation parameter becomes large (which corresponds 
to a very light and/or very thin structure in contact with a high density fluid). During 
this study multiple resonance processes were observed where a single in vacuo resonance 
mode can have several resonance frequencies under heavy fluid loading. The spectrum 
of the operator can no longer be described in terms of the discrete spectrum based on a 
countable series of resonance mode/resonance frequency pairs. This typically nonlinear 
phenomenon is closely linked to the concept of non-linear frequency modes p], where 
the non linearity depends on the time parameters (such as those involved in porous 
materials) rather than on the more usual geometrical parameters. This frequency non- 
linearity is related to problems which are written in the form 0{uj)U = S, where u 
is the unknown factor and S is the source term. 0{uj) is an operator which depends 
non linearly on the angular frequency u, as in the case of porous environments or in 
the context of vibroacoustics where the coupling depends non-linearly on the frequency 
via the Helmholtz equation Green's kernel. When the coupling is weak (as in air), this 
process is not easily observed. 

Here it is established that this behavior is a general property of fluid-loaded plates. 
In paragraph [2], the approximate method for calculating the spectrum of a fluid loaded 
vibrating structure based on a high order perturbation expansion is recalled, and this 
multiple resonance is property is then illustrated by giving a numerical example which 
conflrms the validity of the method (based on comparisons with a numerical resolution 
of the exact equations). Paragraph [3] gives a theoretical description of multiple resonance 
in the case of the simply supported plate loaded with a fluid with an arbitrarily high 
density. It is then established, using classical results on the connection between the order 
of an entire function and the distribution of its zeros, that any in vacuum resonance mode 
in this structure will have an inflnite number of resonances under heavy fluid loading. 
The conclusions and the possible extensions of this work are presented in paragraph |H 



2 HIGH ORDER PERTURBATION METHOD 

Let us flrst recall the main result obtained in our previous study [8] . Consider a thin 
flnite elastic structure (with surface density pph and thickness h) occupying a domain 
S, the behavior of which is described by a differential operator A. This structure is 
loaded with a perfect fluid at rest extending to inflnity. Using perturbation methods, one 
calculates the eigenmodes Um{M, lu), where M is a point of the domain S, and eigenvalues 
Am{uj) = PphCj'^{(jj) as solutions of the classical integrodifferential equation 



where e = pf/pph is a small parameter in the case of a metal structure in contact 
with air (density p/). G{M\ M' ,uj) is the Green's function of the Neumann problem 
in the Helmholtz equation outside the surface S. One deflnes the radiation impedance 
j3^{U^ V) = /s /s U{M')G{M; M', uj)V*{M)dMdM' , where U and V are function deflned 
on the surface S. The weak formulation of the problem then reads flnd Ifm such that for 





every V a{U^, V*) - A^(u;) ((f/^, V*) - ep^{Um, V*)) = 0, where a{U, V*) stands for the 
potential elastic energy of the structure, {U, V*) is the usual inner product, pphuj^{U, V*) 
is the kinetic energy and pphuP'fiujiUm, V*) is the energy lost by radiation into the fluid. 
It is worth noting that because of the fluid-loading, each eigenmode and each eigen- 
value depends on the frequency. Until necessary, this dependence is omitted to facilitate 
reading. 

Let us expand into a Rayleigh-Schrodinger perturbation series [9] both the eigenmodes 
Um and the eigenvalues in e : 

Um{M) = ujfi\M) + eU!^\M) + ... + e^Uj::\M) + ... (2) 
= A(^) + eA« + --- + e^A(^) + --- (3) 

The Rayleigh-Schrodinger Method is one of the most powerful methods available to deal- 
ing with vibroacoustics problems involving low frequency dependent damping such as 
those arising in viscosity or sound radiation in light fluid. Substituting the two perturba- 
tion expansions into the weak formulation given by equation ([T]), collecting and equating 
the coefficients of equal power e to zero yields: 

6°:a(L/i°),\/)-AW(t/i°),K) = 0, (4) 
■■ a{U^\V) - A^2HU!^\V) = A« (f/^, V) - AW/3.(f/i°), \/), (5) 
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Equation (jll) is the usual relation giving the eigenmodes of the elastic structure in 
vacuo. If the eigenmodes f/^-* (M) of the elastic plate in vacuo are normalized, it is easy 
to obtain the first order term of the eigenvalue form equation ([5]), that is 

A(^^ = A(^)/?™™, (7) 

in the previous equation, it has been denoted = /3^(f7^\ t/^^^*). To simplify the 
reading this notation will be used in the rest of the paper. IJ^\M) is calculated as 
an expansion into a series of zeroth-order modes: [/^^(M) = '}2n(^mUl^\M) where the 
coefficients are easy to calculate [9]. One has 

A(o) 

In m nmn /q\ 

"m ~ T(0) a{0)'^'-' ■ 
ivjTt -/in 

For all s > 1, one has: 

A^:) = 'E'A«/5.(t/r'-),^i°)) - e"a«(?7(:-'),^W). (9) 

1=0 1=1 

In a similar manner, each component of the eigenmode's perturbation expansion is de- 
veloped into a series of zeroth-order modes: 

f>i^nM)=E«™^l°nM). (10) 



One shows that the expansions of the eigenvalues are given by: 
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Now let us calculate the eigenvalue expansion up to order 5: 
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Upon re-ordering the terms in equation (|T5|) . one obtains: 



a(o) 



1 + e/?:r" + e'(/?"'")' + e'(/?:r")' + e'(/3:r")' + e'(C)' + 
+e' E [l + 2e/?™™ + 3e2(/?:r™)2 + 4e3(/?™-)3 + 

E [i + '^epr + 362(c"')2 + ■ • 

PT^m 



<! E [1 + ^epr + ■■■]+ E I [i + se/^:^'- + ■ 
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(16) 



For e/^™"* < 1, one can recognize in each square bracket in equation ( 1161) the expansion 
of 1/(1 — 6/5^7"*) its successive powers. By identifying the expansion with the original 
functions, one obtains : 



jm - (1 - 6/5--) ^ ' (1 - 6/5-^)2 + ' (1 - 6/3--)2 
5 



+e^--- (17) 

Now by limiting this expansion to the first term in equation (ITTj) . one obtains the 
high order perturbation expansion: 

A™, 1 



AW (l-e/5--)- ^ ^ 

It is worth noting that while written as a function of e alone, the relation (fTSi) can not 
be considered as a first order expansion since it takes into account an infinite expansion 
into powers of e. 

Now by introducing the pulsations, one obtains the expression of the high order, 
frequency dependent, eigenpulsations: 

~ ^^Tt. — \ — ^- (19) 

™ (1 ) 



It can be shown [8j, that the high order perturbation expansion given in equation ( 1191) 
remains valid even for e ^ 1. 

However, if one is interested in the resonance modes Um{M) and resonances pulsations 
LJm, given by A^ = pphu'^, that are the non-vanishing solutions of equation a{Um, V*) — 
pphw"^ (^{Um,y*) — e/5(im(f^m5 ^*)) = calculated using perturbation expansion (|T9|) . it 
can be easily seen that one has to solve the following resonance equation 

= ^li^m) ^^i^ ^Ti (i - e/^r ) ^ ^m (i - ^cr) ^ ^^T- (20) 

This calculation is extremely difficult because one needs to find the complex roots of a 
non convex complex function. When attempting to solve this equation, the following 
cases occur: 

• For e = 0, there is obviously one real resonance (and its opposite) solution of 

ul = ul = Cj^T- (21) 



• For e -C 1, the fiuid loaded resonance pulsations are very close to the in vacuo ones. 
And it can be shown that there there is one complex resonance (and the opposite 
of its conjugate —Cj'^ solution of the approximate equation in which the modal 
impedance is estimated at the in vacuo eigenpulsation: 

this relation is thereafter named the classical perturbation expansion resonance 
equation. 
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Figure 1: Comparison between the real (left) and imaginary (right) parts, both in Hertz, of the complex 
resonance frequency / of the mode (m = 1, n = 3) in the case of a steel plate in contact with water when 
varying its thickness h in meter, obtained by using the exact solution (thin curves with circles) and the 
high order perturbation expansion (thick curves). 



• For e = 0(1), there are one or more complex resonances (and the opposite of the 
corresponding conjugates), as shown in figure ([1]) that are solutions of 

(^i = co^Ti (i - e/^r ) , 

or, equivalent ly 

(l - e/^r ) - ^'-T = 0, (23) 

this relation is thereafter named the high order perturbation expansion resonance 
equation. 

The corresponding resonance modes are simply given by the calculating the eigen- 
modes at the corresponding resonance frequency Um{M) = Um{(^m,M). 

As an example, in the case of a clamped steel plate measuring Im x 0, 7m in water, 
the resonance frequency in the case of the mode m = 1 n = 3 calculated by solving 
the exact equations using a Boundary Element Method can be compared in figure [1] 
with those obtained using the high order perturbation method given by the numerical 
resolution of the high order perturbation expansion resonance equation fl23l) when varying 
the thickness of the plate. It is worth noting that the exact solution presented in figured] 
required approximately 1 month to calculate on a 4 processor parallel computer whereas, 
the calculation of the approximate solution required only a few minutes on the same 
computer using Mathematica [6]. In a vacuum, the real part of the frequency increases 
with thickness, while the imaginary part is zero; also as lj^-'^ oc h"^ and e oc l/h, the 
function tu^ — uj'^'^ / (1 — ^13^^) ~ oj"^ for /i ^ 0. At a given thickness, the mode m = 
1 n = 3 has at least two resonance frequencies located on either an "upper branch" 
(discontinuous curves) or a "lower branch" (continuous curves). It is worth noting that 
this behavior occurs at all the resonance modes of this plate. 

These results show not only that the behavior of this mode can be described very 
precisely using the high order perturbation expansion, but also (and this is the key point 
in the present paper) that a given resonance mode can have several resonance frequencies. 

3 MULTIPLE RESONANCES 

Here we use the analytical expression for the resonances given by the high order 
perturbation expansion that is the various roots of equation 0231) . that is the complex 



pulsations a; = a;™, such that u'^{l — e/?^) = cD^-*^. The whole problem focus here on the 
existence and the counting of these roots. 

The questions arising about the existence and counting of roots can be answered by 
using some classical results based on the theory of entire functions pUl |3]. There is 
a strong link between the distribution of the zeros of entire functions and their order. 
Let us recall that if f{z) is an entire function, the function M{k) = maxi^i^;, \ f{^)\ 'wiH 
increase indefinitely together with k. The rate of growth of M{k) can be characterized by 
comparing it with function exp(A;). By definition, the order (of growth) p is the smallest 
number such that M{k) < exp k^^^ for every 5 > and for k > ko(S). From this, it can 
be easily seen that the order is given by the formula p = limsup^^(,olnlnM(fc)/lnfc ; 
as examples, a polynomial in z is of order and the exponential function is of order 1. 
Moreover, the Jensen formula ([lOl H]) states that if nf{k) is the number of zeros of an 
entire function f{z) of order p inside a disk with a radius k centered at the origin, then if 
k is sufficiently large, one has nf{k) < kP+\ 6>0 (see [10], theorem (8.3) p.328). Since 
nf{k) is a non- decreasing function of k it has either an infinite number of zeros or none 
(as in the case the function exp(z), which is never zero and is of order 1). 

Now the key point here is how to determine the order of the complex function f{k) = 
k"^ — k"^ — ek'^l3"^{k) = 0, with /3™(A;) = /J™*" and km if a fixed real number. Since k"^ — k"^ 
and k"^ are zero order functions, using theorem (6.4) in [10]: if Fi{z) and F2{z) are entire 
functions of order pi and p2, respectively, and if pi < p2 then the order p of the sum 
Fi{z) + F2{z) is equal to p2 and theorem (6.7) in [10] :if F{z) is an entire function of order 
p, and P{z) is a polynomial of positive degree, then the product F{z)P{z) has order p, 
it is easy to see that the order of f{k) depends only on the order of P"^{k) when e 7^ 0. 

Now the next point is to determine the order of the radiation impedance j3™'{k) as an 
entire function of k. For this purpose, one shows with the simple example of a simply 
supported plate (this case was chosen because of the relative simplicity of the calculations 
involved) that each mode has an infinite number of resonance frequencies. 



3.1 The Simply Supported Plate 

Now let us take the more complex case of a simply supported rectangular plate (di- 
mensions a X b) vibrating in an m x n mode. In a vacuum, the mode is given by 
U^}^{x,y) = 2 / \/ab sin mirx/ a sin niry/b. The key point is again how to determine the 
order of the radiation impedance. The procedure developed here is slightly different 
from the previous one because is seems to be impossible to construct an entire radiation 
impedance series in terms of to the variable ; it is based on a method involving com- 
binations of entire functions [TU]. By making a change of variable it can be easily 
shown that PkiUj^l Uj^l*) = /3'""(fc) is given by 

db ^tkD 
(3'^^{k) = / / / / sin mvrx sin riTT?/ sin mvra;' sin nvr?/' dxdydx'dy' (24) 

TT JO Jo Jo Jo TtD 
Anh /■! pikR 

= -— GUX)G^iY)—-dXdY, (25) 

TT JO JO TTK 



where D = ^a^{x - x'Y + 62(y - y')2 and R = y/a^X^ + bW^. The functions G'^(x) 
and Gn{y) are given by single analytic integrals 

Gmix) = / sinmnix + x') sin mnx'dx' (26) 

JO 



sin miTX + m7r(l — x) cos mirx 
2rmT 

Now let us make a new change of variables from rectangular to polar coordinates 

4 r „ /x\ „ /v\ e*'"" 



(27) 



with 



= — / Jlmn{k,e)de - - J21mnik,e)de, (29) 

TT JO TT Je„ 



J>.„„(»=.<') = /;^^G„(!:^)G„(l^;i^)e-* (30) 



_ fr cos 6* \ _ / r sin ' 



J2mnik,e) = G^l^-^jGni^-^je'^'-dr (31) 

The closed formulas obtained for both Jimn{k,6) and J2mn{k,9) using Mathematica [6] 
analytical integration contains thousands of terms (they are about 850 A4 pages each), 
making it impossible to be controlled or used. Then an approximation of these functions 
that preserves they order property of f3"^^{k) had been build. This approximation is 
denoted by /3'"'^(fc). 

To construct this approximation, one can remark that the function to be studied is 
the sum of functions with of following form: 

fdc rg{B) 

F{k) OC / f{r,9)exp{ikr)drd9, (32) 

JO JO 

The exact expression of g{6) does not matter, what is important is that this function 
does not depend on k. Now by using integration by parts with respect to r one obtains 

F{k) oc f(g(9),9)—^d9 + Oik') (33) 
f(a(9) 9) e*^f(^) 

oc I ^.(jj 9'me + 0{e), (34) 

again by using integration by parts with respect to 9 one obtains: 

pikg{ea) 

F{k) oc H{9^)^^ + 0{e), (35) 

with H{9) = f{g{9),9)/ g'{9). Then the order of the function F{k), which is given 
by the rate of growth of the function M{k) = ma:X\k\^ oo F{k), is given by M{k) = 
meiX\k\-,oo H {9a) exp{zkg{9a))/k'^, then the order does not depend on a exact description 
of the functions H{9a) and g{9a) and the approximations can be made on these functions 
without changing the impedance order. 

Since one always has < 9a < tt/2, one can expand the functions to be integrated 
with respect to the angular variable 9 up to the order 1: Jimn{k,9) around 9 = and 
J2mn{k,9) around 9 = 7r/2. This gives 

^, ,,, ta'^k^rriTT + 2m^7r^ — lakm^TT^ — (—l)"^2e^"'''m^TT^ , , 

^Imn^,^ -Jlmfc = a — 36 

4rmT{ak — rmTj^[ak + rmTj'^ 

J2mn{k,9)^32nik) = b -— ^ 4 37 



ji,2m{k) are regular entire functions of k, in particular, it can easily be shown that 



jim{k = ±imr a) = o s o (38) 

3Uk = ^) = a (39) 

,.„(fc = ±n./6) = ,6,.,3,. (40) 

^-(^ = 0) = & (41) 

Then, by integrating with respect to 9, one obtains P'^"'{k) = —4:/7T{jim{k)9a + 
32n{k){n/2 — 9a))- It is worth noting that the accuracy of the approximation of /9"'"'(/c) 
decreases as the mode order increases. However, at low orders (say m < 3 n < 3) the 
accuracy is very good. Now as jim{k) and j2n(fc) are entire functions of k, /3'""(A;) is also 
an entire function of k. It still remains to estimate the order of /3™"(fc). For given m and 
n and for \k\ oo, one has 

-(-l)"2nW'*= , , 

J2n{k) ^ ^ , (43) 

r-{k) ^ ^ ((-l)"^^«e-'^ + (-1)"^ - 9a) eA 1. (44) 



and 



The order of /?™'"(fc) is easy to obtain from classical results on the sum and product 
of entire functions (see for example |10], chapter VII paragraph 6). First it is worth 
noting that the entire function {—l)'^m?/a?'9aQ^"'^ + (— l)"'n^/6^ (7r/2 — 9a) e**^ cannot be 
identically zero for any value of a, b, m or n when \k\ —>■ oo. To see this more clearly, one 
can deal separately with the cases a = b and a.y^ b. 

The case a = b. One has 9a = 7r/4 and /3™"(A;) ^ ^ {{-l)"'m^ + (-1)"^^) e*°^^ 
becomes the product of the exponential function e*"'^ of order 1 and the inverse of a 
degree 4 polynomial in k, which is of order 0. $"^^{k) is therefore an entire function of 
order 1. 

The case a ^ b. Let us take for example a > b. Then with \k\ oo one has 



max 



-.lak 



^ max 



it is worth noting that this is true in the half plane '^{k) < 0, 
whereas in the half plane 53(fc) > 0, $"^"'{k) tend rapidly to zero. In the half plane 
^{k) < 0, one therefore has /3''""(fc) ^ ^ ((— l)™^6'Q,e*"'^^ ^ which is also the product 
of the exponential function e*"'^ of order 1 and the inverse of a degree 4 polynomial in k, 
which is of order 0. $^"{k) is therefore an entire function of order 1. 

From the Jensen formula, the number of zeros in an entire function f{z) or order 1 
inside the disk with radius r centered at the origin, if r is sufficiently large, is given by 
n{k) < k^^^ , 6 > and since n{k) is a non-decreasing function of k it has either an 
infinite number of zeros or none. Then, with \k\ — > oo, the radiation impedance of the 



Fi gure 2: Contour plots of the amplitude of the radiation impedance in the mode m = 1, rt = 1 in the 
case of a simply supported plate (1 m x 0.7 m ). Left: approximate /3^^{k), right: exact f3^^{k). 

membrane has either an infinite number of zeros or none. To estabhsh that there is at 
least one zero (and then an infinite number), one apphes the argument principle whereby 
if C is a disk with radius K centered on the origin, the number of zeros inside this disk 
n^{K) is given by the integral np{K) = 1/(2ot) //3ik)dk, with k = Ke'^ . 

The number of zeros inside this disk, in the case of the mode mn, denoted by n^mn (-ft'), 
is given by the integral n{K) = l/(2m) §c $'"''%k) / ^^'^'\k)dk, with k = Ke'\ if /3"'"(A;) 
has no pole inside C. For example, with a = 1 and 6 = 0.7, one computes the previous 
integral numerically. It can easily be established that in the mode m = 1, n = 1 
n^ii(ll) = 0, n^ii(12) = 2, n^ii(18) = 4 and in the mode m = 1, n = 3 n^i3(7) = 0, 
n^i3(8) = 2, n^i3(19) = 4. Again, the modal impedance has an infinite number of 
zeros. Therefore, while in a vacuum, each mode of a simply supported plate has only one 
resonance frequency (and its opposite), under heavy loading condition each mode will 
have an infinite number of resonance frequencies. 

Figures [2] and [3] give contour plots for the amplitude of the radiation impedance 
f3"^^{k) and its approximation /3'""(/c) in the mode m = 1 n = 1 in figure [2] and in the 
mode m = 1 n = 3 in figure [31 As previously, on these curves, the dark spots correspond 
to the zeros of the radiation impedances. It is worth noting that the lower the mode 
order, the better the approximation /3'""(fc) becomes. 

It is worth noting that all the figures [2] to [6] are plotted in a quarter of the complex 
plane. The upper half complex plane has no root: the frequency dependence chosen 
exp{—zu!t) implies a complex pulsation u with negative imaginary part to ensure a de- 
crease of the solution for t — > oo. The curves are plotted on the quarter lower plane with 
either positive or negative real part without any particular difference since the various 
resonances are symmetric with respect to the imaginary axis: resonances occur alway by 
conjugate pairs of the form Um and —uj^ when Um G (D or Um and —Um when G IR 

3.2 Comments About the General Case 

Now let us return to the initial question of the existence or non-existence of multiple 
resonances in fluid-loaded plates. Let us recall that in the case of heavy loading, one 
has to look for the resonance pulsations Umn such that uj^ni^ ~ ^PE'^ni^)) ~ ^mn = 0- 



Fi gure 3: Contour plots of the amplitude of the radiation impedance in the mode m = 1, n — S modes 
in the case of a simply supported plate (1 m x 0.7 m ). Left: approximate /3^^{k), right: exact /3^^(fc). 

As shown in the previous paragraphs, this behavior occurs in a large class of problems. 
Although this statement has not been proved to be be true under all possible boundary 
conditions, it seems reasonable to assume that under boundary conditions of all kinds 
and under very heavy loading condition every resonance mode in a plate has an infinite 
number of resonance frequencies. 

But in the case of light fluid loading conditions, as shown with the previous results 
whatever the loading conditions, the resonance value relation has several zeros. One close 
to the real axis with a small imaginary part a a set of resonance with very large imaginary 
part (see left curve in figure H]). The later are very difficult to detect because of their large 
imaginary part tha implies a very rapid vanishing of their contribution. But as shown 
by Sanchez [11], under light loading condition (as in the case of a steel plate in contact 
with air, for example), there exist only the resonances of the structure alone (one in each 
mode shifted toward the complex plane) and the complex scattering resonances of the 
fluid in contact with the rigid solid. This is an apparent contradictory situation because 
Sanchez deals with a slightly different problem since he looks for scattering frequencies 
for a rigid structure in vacuum with zero eigenvalue (that is cj^^^ = 0) which has only a 
set of scattering resonance frequencies. 

The results presented in left curve in figure H] also show that except for the zero close 
to origin, the zeros of the resonance value relation corresponding to the plate loaded 
with water are very similar to those given by the radiation impedance. In that case, 
contrarily to the light fluid loading situation, both structure resonance and scattering 
resonances have their imaginary part with comparable magnitude, that allows a possible 
identification of the corresponding modes. 

Now under moderate loading condition (as in the case of a steel plate in contact with 
water), what is the relevant relation ? This seems to be an open question. Although the 
resonance relation fc^(l — e/3™'"(fc)) — k^^}^ = is true under very heavy loading, it has to 
be supplemented under small or moderate loadings conditions (that is, when e = C(l)) 
by terms accounting for the cross modal impedance, which make it almost impossible to 
conduct the above analysis. 
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Figure 4: Contour plots of the amplitude of the resonance value equation — e/3™"(fc)) — kmn in 
the mode m = 1, n = 3 in the case of a thin simply supported plate (1 m x 0.7 m x 2 mm). Left: steel 
plate in contact with air (e = 0.08), right: steel plate in contact with water (e = 64). 



3.3 Numerical Example on a Clamped Steel Plate in Water 

To conclude this analysis, one presents in figure [5] and figure [6] a comparison of the 
location of the resonances for a clamped steel plate in contact with water, using contour 
plots. In these curves the various dark spots show a rough estimate of the roots. In 
figure O it is showed a contour plot of the Logarithm of the determinant amplitude 
of the linear system of simultaneous equation obtained by solving the exact Boundary 
Integral Equation Method using a Tchebycheff-coUocation method [7] and in figure [6l 
it is plotted the contour plot of F{lj) defined as a product of the thirty firsts high 
order resonances equations (corresponding to the modes (1,1) to (6,5)) for the same 
plate: F{uj) = n^^=\ Il^^^e^^ (1 - ^PT) - ^in}? which roots are obviously those of the 
corresponding high order resonances equations. While the absolute values (and contour 
lines) can not be directly compared, the approximate location of the roots of the function 
F{uj) gives a reasonable image of the true resonances. These figures show also that the 
spectrum of this fluid-loaded structure is not a combination of the resonances of the 
structure alone (with frequency lowered by added mass and shifted through the complex 
plane by radiation damping) and of the scattering resonances of the rigid body alone. 
Each scattering resonances set is shifted through the complex plane differently by each 
resonance mode of the structure allowing a very complex distribution of the resonances 
of the fluid-loaded structure. It is worth noting that this distribution of the resonances 
makes a direct numerical estimate of the spectrum quite impossible. This also shows that 
there is always the need of analytical estimates such as high order perturbation expansion 
to be able to construct a numerically efficient estimate of the true spectrum. 

4 CONCLUSIONS 

Here we have shown that the loading of a rectangular plate with a high density fluid 
transforms each resonance into an infinite number of resonances. These results open 
up new prospects. In particular, it would be interesting to obtain similar estimates on 
the vibratory modes occurring in clamped plates or mechanical structures of other kinds, 
although these may not be easy to obtain. Another interesting aspect concerns the asymp- 
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Figure 5: Contour plots of the Logarithm of the determinant amphtude of the hnear system of simul- 
taneous equation obtained by solving the Boundary Integral Equation Method for a clamped steel plate 
in water. 
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Figure 6: Contour plots of the Logarithm of the product of the thirty firsts high order resonances 
equations (corresponding to the modes (1, 1) to (6,5)) for a clamped steel plate in water. 



totic properties of the zeros of entire functions. As we have seen above, making spectral 
estimates is equivalent to finding the roots of the function f{k) = k'^(l — eP"^"{k)) — ^ 
but f{k) is an order 1 function at e 7^ 0, and it can easily be shown that it is also a func- 
tion of normal type [1], then by definition f{k) is an entire function of exponential type. 
Entire functions of this kind which satisfy the inequality /I^j^ln"*" + k'^)dk < 00 

belongs to class A ^ and have regularity properties as regards the asymptotic distri- 
bution of their zeros. Another question which arises is how to determine under what 
conditions moderately coupled systems (such as a steel plate in contact with water) will 
have these properties. Another important aspect worth investigating is the validity of 
modal expansions in the case of resonance modes having two or more resonances frequen- 
cies. It would also be interesting to develop experimental device which would make it 
possible to test the occurrence of the behavior described here. 
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